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(57) ABSTRACT 

Systems and methods for sensing air outside a moving 
aircraft are presented. In one embodiment, a system includes 
a laser for generating laser energy. The system also includes 
one or more transceivers for projecting the laser energy as 
laser radiation to the air. Subsequently, each transceiver 
receives laser energy as it is backscattered from the air. A 
computer processes signals from the transceivers to distin- 
guish molecular scattered laser radiation from aerosol scat- 
tered laser radiation and determines one or more air param- 
eters based on the scattered laser radiation. Such air 
parameters may include air speed, air pressure, air tempera- 
ture and aircraft orientation angle, such as yaw, angle of 
attack and sideslip. 
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OPTICAL AIR DATA SYSTEMS AND 
METHODS 

RELATED APPLICATIONS 

This application claims benefit of priority to U.S. Provi- 
sional Patent Application No. 60/400,462, filed Aug. 2, 
2000, which is hereby incorporated by reference. 

U.S. GOVERNMENT RIGHTS 

This invention was made in part with the support of the 
U.S. Government; the U.S. Government has certain rights in 
this invention as provided for by the terms of Grant #NAS4- 
0243 awarded by the NASA, Dryden Flight Research Cen- 
ter. 

BACKGROUND 

An Air Data System (“ADS”) provides sensed telemetry 
informing pilots, navigators or Vehicle Management System 
computers of air parameter(s) affecting aircraft stability. 
These air parameters include, for example, air speed, air 
temperature and air pressure, each being useful for naviga- 
tion and flight control. The ADS exists in many forms, for 
example, as mechanical, opto-mechanical or opto-electronic 
devices. 

One mechanical ADS includes a Pitot tube, which affixes 
to the exterior of the aircraft. The Pitot tube is a pneumatic 
measuring instrument that receives air external to the aircraft 
and determines air speed based upon air pressure. Atypical 
Pitot tube uses a pressure transducer to measure the pressure 
of air received through a plurality of holes. Determined air 
pressure is combined with air density measurements and air 
temperature measurements using the well-known Bernoulli 
equation to calculate air speed. Pilots and Vehicle Manage- 
ment Systems use this information for decisions regarding 
aircraft stability, flight control and navigation. 

The prior art ADS has limitations. With respect to the 
Pitot tube, for example, if aircraft velocity is too low, the 
pressure transducer may lack requisite sensitivity to provide 
meaningful telemetry data, such that transducer signal noise 
or error may be greater than meaningful differences in 
pressure measurement. Alternatively, if the velocity is very 
high (e.g., supersonic), certain assumptions, such as those 
regarding incompressibility of air flow in Bernoulli’s 
equation, are violated and the measurement is incorrect. 
Accordingly, Bernoulli’s equation must be altered to com- 
pensate for various breached assumptions. This compensa- 
tion process must continue in rapidly changing conditions of 
altitude and air density. It is, therefore, desirable to improve 
accuracy and operability of the ADS. 

One improvement is an Optical Air Data System 
(“OADS”), which uses light to determine parameters of air 
speed. The OADS transmits light pulses into the atmosphere 
and receives light that is reflected, or "backscattered”, from 
aerosols toward the aircraft. Aerosols are fine solids and/or 
liquid particles suspended in air or other cases. 

The OADS may also measure the Doppler effect, receiv- 
ing the backscattered light and measuring return frequency 
to determine speed. However, the prior art OADS relies on 
scattered light that is unpredictable because of varying 
aerosol distributions. For example, aerosol distribution var- 
ies significantly with altitude and cloud content. In addition, 
some regions of the atmosphere contain too few aerosols to 
enable reliable air data measurements. Thus, aerosol-based 
optical air data systems cannot determine the air speed at all 
altitudes frequented by modern aircraft. Finally, an OADS 
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that relies solely on aerosol scattering cannot determine air 
temperature or air pressure. Air temperature and pressure are 
critical air data parameters for determining air density and 
Mach number. Hence, there is a need to improve optical air 
5 data systems and methods. 

SUMMARY 

In one embodiment, an Optical Air Data System 
(“OADS”) advances the art and overcomes the problems 
10 mentioned above by determining air parameters (air speed, 
air temperature and air pressure) based upon laser energy 
that is backscattered from both air molecules and aerosols. 
In one example, the OADS projects laser energy to the air 
and receives backscattered laser energy from the air mol- 
15 ecules and/or aerosols. The OADS converts such backscat- 
tered laser energy to electronic signals which are processed 
to determine desired air parameters. Such air parameters 
may include, for example, air speed, air temperature and air 
pressure, among others, and provide telemetry for air vehicle 
20 flight control and navigation. 

In one embodiment, the OADS includes at least one laser 
configured for generating laser energy. The OADS also 
includes one or more transceivers configured for projecting 
,,, the laser energy to the air. Each transceiver subsequently 
receives backscattered laser energy from the air. Aportion of 
the backscattered laser energy received by each transceiver 
is filtered prior to electronic signal conversion. 

In one embodiment, each transceiver includes a vapor 
50 filter that filters the backscattered laser energy within a 
predetermined band. In one example, a transceiver includes 
a fixed frequency vapor filter, such as an atomic vapor filter 
or a mercury vapor filter, that provides “notch” or “band 
stop” filtering of the backscattered laser energy. The band 
5 5 stop filter may provide peak absorption of background light 
within the band stop region of the filter. In one embodiment, 
the laser is tunable and generates the laser energy at a 
center-tuned wavelength that substantially corresponds to a 
center wavelength of the band stop of the vapor filter. In 
40 another embodiment, the center wavelength of the filter 
corresponds to an absorption region for ultraviolet wave- 
lengths occurring in ozone (0 3 ). Once filtered, a transceiver 
converts the backscattered laser energy into a computer 
processable electronic signal for determination of the 
4 5 desired air parameters. 

In one embodiment, each transceiver includes one or 
more detectors or transducers configured to receive filtered 
and unfiltered backscattered laser energy for conversion to 
electronic signals. In one example, a transceiver includes an 
50 optical transducer that converts a particular range of wave- 
lengths to computer processable signals used to determine 
the desired air parameters. Accordingly, in one embodiment, 
the OADS includes a computer for processing the electronic 
signals to determine the desired air parameters. Other fea- 
55 tures are apparent within the description that follows. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows one Optical Air Data System (“OADS”); 

FIG. 2 shows one OADS; 

FIG. 3 shows one graph useful in illustrating exemplary 
an air speed calculation with an OADS; 

FIGS. 4-9 show graphs illustrating exemplary calcula- 
tions for other air parameters with an OADS; 

65 FIG. 10 shows one OADS transceiver; and 

FIG. 11 shows a flowchart of one exemplary methodical 
operation of an OADS. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows one Optical Air Data System (“OADS") 101 
mounted on or within an aircraft 102. In this embodiment, 
OADS 101 is configured for projecting laser energy 103 
(e.g., laser radiation) to air 104. Impingement of laser energy 

103 on air 104 and aerosol particles 105 (in air 104) causes 
scattering of laser energy 103, which is represented in FIG. 

1 as an energy scatter field 106. Distance between aircraft 

102 and energy scatter field 106 is controlled by the amount 
of overlap between the out-going laser energy 103 and the 
transceiver 110 field of view at a distance from aircraft 102, 
to provide an optimized intensity for return laser energy 107 
and to eliminate the overwhelming measurement error aris- 
ing from displaced air proximate aircraft 102. OADS 101 
detects return laser energy 107 that is backscattered from air 

104 at laser energy scatter field 106. Laser energy 107 may 
be in the ultra-violet (UV) spectrum, for example, having a 
wavelength within a range of 250 nm to 270 nm. Other 
ranges operable to produce scatter field 106 may be used as 
matter of design choice. 

Return laser energy 107 typically contains molecular 
scattered (e.g., Rayleigh) components 107A and/or aerosol 
scattered (Mie) components 107B. OADS 101 distinguishes 
the molecular scattered components 107A from the aerosol 
scattered components 107B and correspondingly determines 
one or more air parameters based on the scattered laser 
energy 107. Examples of such air parameters include air 
speed, air pressure, air temperature and aircraft orientation 
angles relative to the local wind. OADS 101 may be 
configured with other aircraft as well, such as unmanned air 
vehicles (UAVs), helicopters, gliders and space shuttles. 
Although illustrated within a “nose” 108 of aircraft 102, 
OADS 101 may be configured in any other areas of aircraft 
102 . 

As shown in FIG. 1, OADS 101 includes a laser 109 
configured for generating laser energy 103. Transceiver 110 
is configured for transmitting laser energy 103, from laser 

109 via optical coupling 111, and receiving backscattered 
laser energy 107. Optical coupling 111 may exist in the form 
of a fiber optic connection or free space transmission. 
Accordingly, transceiver 110 projects the laser energy as 
laser radiation 103 to air 104. Air 104 scatters laser radiation 

103 at energy scatter field 106 in a plurality of directions 
(e.g., illustrated as vectors 112). Scatter field 106 also 
returns, or backscatters, laser radiation 107 towards trans- 
ceiver 110, which subsequently receives the backscattered 
laser radiation 107. Transceiver 110 converts backscattered 
laser radiation 107 to processable electronic signals, via 
computer 113, to determine the air parameters. 

Computer 113 communicatively couples with transceiver 

110 and processes signals from transceiver 110 to distin- 
guish molecular scattered component 107A from aerosol 
scattered component 107B. Computer 113 determines the air 
parameters based on laser radiation 107 backscattered from 
molecules and/or aerosols in air 104. Accordingly, as 
described below, computer 113 employs one or more digital 
signal processing algorithms to determine such parameters. 

While OADS 101 illustrates one transceiver 110 in an 
exemplary embodiment, a plurality of transceivers may be 
used depending on application. For example, a helicopter 
employing OADS 101 may use two transceivers 110 to 
determine air parameters such as forward velocity (e.g., air 
speed) and a horizontal plane, or “yaw”, of the helicopter. 

An airplane may use three transceivers 110 positioned in a 65 
particular manner to determine various aircraft geometries, 
such as angle of attack and sideslip, in addition to the air 
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parameters of air speed, air pressure and air temperature. In 
addition, air vehicles (fixed wing and rotary) may employ 
three or more transceivers and/or lasers to increase Air Data 
System reliability through a redundant system architecture. 
5 Using three OADS transceivers mounted orthogonally to 
one another may fully resolve a total airspeed vector by 
providing three independent measurements for the air speed 
vector (i.e., corresponding to three axes of a Cartesian 
coordinate system). The transceivers are located in uncom- 
mon planes and their geometry respective of an aircraft 
center-line known. Vector algebra may then be used to 
determine the full airspeed vector, including forward air 
speed, angle -of-sideslip and angle-of-attack. 

FIG. 2 shows one OADS 140. OADS 140 illustrates 
15 another embodiment used for determining air parameters, 
such as those described in FIG. 1, based upon laser energy 
backscattered from both air molecules and aerosols. In this 
embodiment, OADS 140 includes laser 141 configured for 
generating laser energy 142. Laser 141 may be tunable laser 
,, 0 having a tuned center wavelength of about 253.7 nm, 
although other wavelengths may be used. For example, laser 
141 may be a frequency quadrupled, Nd:YAG (i.e., 
neodymium:yttrium-aluminum-garnet) pumped Ti:Sapphire 
(titanium-sapphire) laser. Alternatively, frequency- 
25 quadrupled Yb-doped (ytterbium -doped) fiber lasers can be 
used that offer important benefits of smaller size, lighter 
weight, increased robustness and improved reliability, as 
compared to Nd:YAG-pumped Ti:sapphire lasers. Laser 141 
generates laser energy tunable across 40 GHz in about 100 
, 0 MHz spectral steps. Additionally, laser 141 may tune +/-20 
GHz about a center frequency of approximately 1182.5 THz, 
or c/253.7 nm, where c is the speed of light (approximately 
3x10 s m/s). Laser 141 radiates laser energy 142 to beam 
splitter 143, which splits the beam into two components 
, 5 143A/143B. Component 143Ais directed through air 144; 
component 143B is directed to beam splitter 145. 

In particular, component 143A of laser energy 142 
directed to air 144 is scattered into scatter field 146. Scat- 
tering of component 143Ais illustrated by scattering vectors 
40 147 in scatter field 146, whereas return scattering is illus- 
trated by backscattered laser radiation 148. Component 
143B of the laser energy 142 is used as a reference for 
comparison to backscattered laser radiation 148. Such a 
comparison is for example useful in determining air param- 
45 eters such as air speed, since transmitted and received 
frequencies of the laser energy can be ascertained for use in 
a Doppler equation; such a process is explained in greater 
detail herein below. 

In the illustrated embodiment, backscattered laser radia- 
50 tion 148 is received through optics 149. In one example, 
optics 149 is a telescope that demagnifies backscattered laser 
radiation 148 into a laser beam 150. Optics 149 also directs 
beam 150 to beam splitter 151, to split beam 150 into two 
components 150A/150B. Component 150B of beam 150 
55 passes through vapor filter 152 to detector 153 to produce 
electronic signal 158 representative of the component 150B 
impinging detector 153; whereas component 150A is 
directed by beam splitter 151 to detector 154. 

In one embodiment, detector 154 is a photodiode that 
60 receives laser 148 radiation and converts such radiation into 
an electronic signal 155. Detector 154 connects to a central 
computer 156 to process electronic signal 155. Similarly, 
detector 153 is a photodiode configured for detecting com- 
ponent 150B, which is filtered by vapor filter 152 as filtered 
component 157. Detector 153 converts component 157 to an 
electronic signal 158 for processing by central computer 
156. 
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Accordingly, electronic signal 158 corresponds to back- 
scattered laser energy 148 as filtered by vapor filter 152; and 
electronic signal 155 corresponds to unfiltered backscattered 
laser energy 148. Electronic signal 155 is thus used to nullify 
certain anomalies as computer 156 processes electronic 
signal 158. For example, when processed with electronic 
signal 158, signal 155 may be used to remove, from signal 
158, certain laser transmission power fluctuations in filtered 
component 157 caused by atmospheric changes in air 144. 
Such a process is explained in more detail in FIGS. 4-7. 

Reference component 143B of the laser energy 142 is split 
into two components 159/160 by beam splitter 145. Com- 
ponent 160 is directed by beam splitter 145 to vapor filter 
152 via mirrored surface 161, to measure filter 
characteristics, whereas component 159 is directed by beam 
splitter 145 to detector 162, to generate electronic signal 
163. Electronic signal 163 is for example used to normalize 
power fluctuations in the return of backscattered laser radia- 
tion 148 caused by power fluctuations in the generation of 
laser energy 142 by laser 141. Such a process is explained 
in more detail in FIGS. 4-7. 

Vapor filter 152 filters component 160 to produce filtered 
component 164. Filtered component 164 is directed to 
detector 165, via mirrored surface 166, and then converted 
to an electronic signal 167. Central computer 156 processes 
electronic signal 167 to determine filter characteristics, such 
as frequencies and suppression features of the band stop 
region of vapor filter 152. Such a process is also explained 
in more detail in FIGS. 4-7. 

It should be noted that FIG. 2 displays the OADS 140 
using free space optical transmission and optical 
components, such as beam splitters 143, 145 and 151, and 
mirrors 161 and 166. However, a more environmentally 
robust method may use optical fiber for laser 141 energy 
transmission along paths 142, 143A, 143B, 159, 160, 164, 
150, 150A, 150B and 157. In such an embodiment, fiber 
splitters are used for 134, 151 and 145 and mirrors 161 and 
166 are eliminated. 

FIG. 3 shows one graph 200 useful in illustrating an 
exemplary air speed calculation with OADS 140. Graph 200 
shows two curves 201/202 comparing normalized layer 
energy magnitudes as a function of frequency (normalized 
laser energy magnitude plotted with respect to axis 205 and 
frequency plotted with respect to axis 204). Curve 202 
exemplifies filtered radiated laser energy such as that of 
filtered component 164 of FIG. 2. As such, curve 202 shows 
filter characteristics of vapor filter 152 of FIG. 2 determined 
by processing of electronic signal 167. Curve 202 shoes a 
peak absorption of filter 152 occurring at a down-translated 
frequency of 0 GHz. By way of example, the actual peak 
absorption frequency of filter 152 may be about 1182.5 THz 
(i.e., having a corresponding wavelength of about 253.7 
nm). 

Laser energy 142 generated by laser 141 passes through 
filter 152 to provide filtered component 164. Once filtered 
component 164 is converted to electronic signal 167 by 
detector 165, computer 156 analyzes and stores features of 
vapor filter 152 through digital signal processing of signal 
172 (e.g., reference features under controlled conditions for 
use in future calculations). As shown in this example, 
features of vapor filter 152 have approximately 10% nor- 
malized absorption at approximately +/-5 GHz (i.e., 0.9 
normalized transmission factor at approximately +/-5 GHz 
according to axis 205) about the peak absorption frequency. 
Other types of suitable filters may include different 
absorption/transmission features. 
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Curve 201 exemplifies filtered backscattered laser energy 
such as that of filtered component 157 of FIG. 2. In one 
embodiment, curve 201 is used to determine air speed by 
comparison to curve 202. For example, curve 202 illustrates 
5 how vapor filter 152 affects laser energy 142; curve 201 
similarly illustrates how vapor filter 152 affects laser energy 
142 as laser energy 142 is backscattered (e.g., returned laser 
radiation 148) from air 144. Frequency shift 203 represents 
the change in frequency of peak absorption for vapor filter 
to 152 between transmitted laser energy 142 and returned laser 
energy 148. Computer 156 processes algorithms applying 
Doppler velocity equation to determine air speed from 
frequency shift 203. 

To determine air speed in one embodiment, computer 156 
i5 determines how far in frequency the peak absorption fre- 
quency of filtered component 157 has shifted from the initial 
laser frequency by comparing curve 202 to curve 201 (e.g., 
comparing peak absorption frequencies of filtered compo- 
nents 157/164). The frequency shift 203 between the two 
20 frequencies substantially equates to a radial wind velocity 
through the Doppler velocity equation: 
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where Av D represents the Doppler frequency shift, v R rep- 
resents velocity component of the vehicle (e.g., aircraft 101 
of FIG. 1) along the laser direction of propagation 143B and 
50 X represents the wavelength of laser energy 142. 

In one embodiment, the wind velocity can be measured by 
determining the frequency shift from curve 202 of graph 200 
as compared to curve 201 of graph 200. This is accom- 
plished by calculating the symmetry point of both spectra 
, 5 202 and 201 and determining the difference in symmetry 
points between the two curves. 

FIGS. 4-9 show graphs illustrating exemplary calcula- 
tions for other air parameters with OADS 140. For example, 
after determining frequency shift due to air speed as shown 
40 in FIG. 3, other air parameters, such as air temperature and 
air pressure may be calculated. In one example, computer 
156 initially determines an intensity measurement of the 
detected backscattered laser radiation (e.g., filtered compo- 
nent 157 detected by detector 153) from electronic signal 
45 158. This experimentally verified intensity measurement of 
returned laser energy corresponds to the following equation: 

Ss{'')=PJ L D s T R E s \dv r \dvl i Uy l )F(v r -v)(rR(v r -{v r \v D ))+ 

mM(v- Av d )))] (Eq.2) 

50 where S s (v) is electronic signal 158 from detector 153; is 
the laser power, T^ is the transmission coefficient through air 
144 along laser path 143B, L(v,) is the laser line shape 
inherent to the laser 141 output as a function of laser 
frequency v ; , T s is the transmission coefficient through air 
55 144 along laser path 148, E 5 is optical efficiency of the 
detector channel through detector 153, F(v) is the band stop 
frequency range of vapor filter 152 centered at a frequency 
of v, R is Rayleigh scattering as a function of frequency 
(applicable to the Rayleigh regime) v r for backscattered 
60 laser radiation minus the quantity of laser frequency v, 
minus the Doppler shift Av D , r is the Rayleigh scattering 
magnitude coefficient dependent on air density and the 
Rayleigh backscattering coefficient, M is Mie scattering as a 
function of v ; . minus the quantity of v ? minus Av^,, m is the 
65 Mie scattering magnitude coefficient dependent on aerosol 
concentration and the Mie backscattering coefficient, and D s 
is detector 153 efficiency of the detector channel 158. The 
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Rayleigh backscattering coefficient r and the Mie backscat- 
tering coefficient m are constant for a particular atmosphere. 
These coefficients correspond to the number of scatterers 
(i.e., molecules for Rayleigh, aerosols for Mie) per unit 
volume of atmosphere. 

With a result obtained for the measured intensity of the 
returned laser energy, computer 156 commences the process 
of determining the other air parameters. Such a process, for 
example, begins with determining characteristics of vapor 
filter 152 via the transmission of reference laser energy 
through vapor filter 152 (e.g., component 160). For example, 
measuring the band stop characteristics of vapor filter 152 
with laser 141, (e.g., via component 143B to electronic 
signal 167) during experimentation yields a convolution of 
the laser wavelength and the filter according to the following 
equation: 

S f (v)=P l E f D f \ dvi_L(v,)F(v r v) 1 (Eq.3) 


In one embodiment, the absorption/transmission charac- 
teristics of the vapor filter 152 are normalized using Eqs. 3 
and 4. Eq. 3 yields the stop band characteristics of filter 152 
and Eq. 4 accounts for power fluctuations in the generation 
of laser energy 142. With the power fluctuations of Eq. 4 
substantially removed, a "normalization channel” is 
revealed and power fluctuations attributable to atmospheric 
changes can be accounted for. 

In one embodiment, additional power fluctuations caused 
by atmospheric changes in air 144 are also removed. For 
example, laser energy detected by detector 154 (e.g., via 
component 150A) assists in removing laser power fluctua- 
tions caused by atmospheric changes in air 144. 
Accordingly, detector 154 converts received laser energy 
into electronic signal 155. Computer 156, in turn, processes 
electronic signal 155 to determine the normalized laser 
energy magnitude according to following equation: 


where S F (v) is signal 167 from detector 165 as a function of 
frequency v (e.g., as illustrated in curve 122); E F is the 20 
optical efficiency of filter 152 collection along paths 160 and 
164, and D F is the optical efficiency of detector 165. 

Note that all optical efficiencies E F and E^ capture signal 
losses that are optical in nature. For example, the optical 
efficiency for detector 165, D F , includes the optical beam 25 
splitting ratios for beam splitters 143 and 145, the optical 
transmission and coupling across filter 152 and the optical 
delivery efficiency onto detector 165. The optical collection 
efficiency for detector 153, E^, includes the collection effi- 
ciency of telescope 149, the optical coupling efficiency into .to 
path 150, the beam splitter ratio of beam splitter 151, the 
transmission efficiency across filter 152 and the delivery 
efficiency onto detector 153. Detector efficiencies D F and D 5 
include the detector conversion efficiencies for detectors 165 
and 153, respectively. Thus, D F is the conversion efficiency 35 
whereby detector 165 converts laser energy along path 164 
into an electrical signal 167. Likewise, D 5 is the conversion 
efficiency whereby detector 153 converts laser energy along 
path 157 into an electrical signal 158. 

Backscattered laser radiation 148 can include power flue- 40 
tuations that are caused by laser 141 while generating laser 
energy 142. Accordingly, laser energy detected by detector 
162 (e.g., via component 159) assists in normalizing power 
fluctuations attributable to laser 141 by removing such 
power fluctuations from the analysis. In one embodiment, 45 
detector 162 converts the laser energy into electronic signal 
163. In turn, computer 156 processes and normalizes accord- 
ing to the following equation: 

S L =P L E L D L jdvL(v), (Eq.4) so 

where S L is the electronic signal 163 from detector 162, E L 
is the optical collection efficiency for detector 162, D F is the 
conversion efficiency of detector 162 and P F is the power of 
laser 141. Note that the optical collection efficiency E L 
includes the beam splitting ratios of beam splitters 143 and 55 


S N =P L T L T a E N D N \ dvj dv{L(v,)(rR(v-(v r Av D ))+ 

mM{v-(v r Av d )))], (Eq.5) 

where is the signal 155 from detector 154; E^ is optical 
collection efficiency of the detector 154 and is the 
conversion efficiency of detector 154. 

In one embodiment, it is advantageous to normalize the 
various characteristic functions to enable a closed-loop 
solution to the process of determining temperature and 
pressure. In one example, therefore, computer 156 calculates 
the normalized laser line shape according to following 
equation: 

[L(y,)dv r l, (Eq.6) 

where (as before) is laser line shape frequency and L 
denotes the laser line shape as a function of frequency. In 
another example, computer 156 calculates normalized Ray- 
leigh Function according to the following equation: 

lR(v r )dv=l, (Eq.7) 

where R denotes the Rayleigh line shape as a function of 
frequency v r , applicable to the Rayleigh regime. In another 
example, computer 156 scales the electronic signal 167 
recorded from detector 165 by dividing all recorded values 
by the maximum value according to the following equation: 

MAX(S f (v))« 1, (Eq.S) 

where MAX denotes an operation that finds a maximum 
value of a particular function, and S F denotes electronic 
signal 167 measured from detector 165, as a function of 
frequency v (e.g. as illustrated in curve 221 of FIG. 4). In 
another example, computer 156 normalizes the Mie Func- 
tion according to the following equation: 

M(v)=d(v), (Eq.9) 


145 and the delivery efficiency of laser beam path 159 onto 
detector 162. 

Curve 221 of FIG. 4 represents the magnitude of laser 
energy (component 164) filtered by vapor filter 152 and 
normalized between 0 and 1. Curve 221 represents the 60 
magnitude of the laser energy as a function of frequency 
(i.e., laser energy magnitude plotted with respect to axis 222 
and frequency plotted with respect to on axis 223). Curve 
221 , therefore, illustrates filtered laser energy via component 
160 as determined by computer processing of electronic 65 
signal 167, plotted as laser energy magnitude normalized 


where 8(v) is the delta function. 

In one embodiment, dividing the signal 167 collected 
from detector 165 (and represented by Eq. 3, above) by the 
signal 163 collected from detector 163 (and represented by 
Eq. 4, above) removes laser 141 power fluctuations, as 
follows: 

Sf(v) PlEfDf J ' dv lW v l) F ( v l~ v )\ <Eq ’ 101 

5ZN) " P L E L D L fdvUv ) 


between 0 and 1 versus frequency. 
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Equation 10 can be simplified to: 

s>(v) EfDf (Eq. 11) 


In one embodiment, tuning the laser 141 to a frequency 
(v A ) far enough removed from the effects of the vapor filter 
152 enables the measurement of the ratio of the optical and 
detector efficiencies of the signal channels 167 (S^, repre- 
sented by Eq. 3 above) and 163 ( S L , represented by Eq. 4 
above). This, in turn, enables the normalization of the signal 
167 measurement to one, for simultaneously checking for 
laser, detector and filter abnormalities on a scan-by-scan 
basis: 

S F (v„) _ EfDf (Eq. 12) 

■sIoT) ~ ~ElDl 

In one embodiment, LF(v) are determined to generate a 
look up table of the convolution of theoretical Rayleigh 
functions (calculated in terms of temperature and pressure) 
with the measured filter function. Since, the measured filter 
function is already the convolution of the laser and filter 
spectra, convolving the Rayleigh function with the measured 
filter signal 167 yields the expected return signal from an 
atmosphere of pure Rayleigh scatterers. 

In one embodiment, the measured signal 158, which is the 
backscatter return from the atmosphere 144, that passes 
through the vapor filter 152 (and is represented by Eq. 2 
above), is divided by the signal 155, which is the backscatter 
return from the atmosphere 144, that does not pass through 
vapor filter 152 (and is represented by Eq. 5 above). This 
removes changes in signal transmission that are independent 
of the factors to be measured: 


10 

and pressure from the return signal. Initially, computer 156 
uses theoretical Rayleigh functions which are functions of 
temperature and pressure in conjunction with the measured 
filter transmission to generate a lookup 170 table of laser, 
5 Rayleigh, and filter (LRF(v)) convolutions that are depen- 
dent on atmospheric temperature and pressure. Computer 
156 then compares the normalized return signal to the 
lookup table 170 to determine atmospheric temperature and 
10 pressure. In order to compare the return signal with the 
lookup table 170 computer 156 accounts for the magnitude 
of the Mie scatterers as well as any changes in air density, 
which changes the magnitude of the Rayleigh signal. 

The vapor filter is used as a bandstop filter for frequency 
I s stability, optical depth, and optimal filter shape. For the 
purposes of separating the Rayleigh and Mie scattering, an 
optical depth of approximately 60 dB provides excellent 
absorption of Mie scattering within a small frequency vari- 
ance around v 0 (i.e., where w-is a normalized frequency of 
20 0 GHz). For example, the 60 dB of absorption in the atomic 
vapor cell exists in a region that is not contaminated by Mie 
scattering. This 60 dB absorption region is used in acquiring 
initial estimates of pressure and temperature (explained 
below in FIG. 5). Such absorption is observable below in 
FIG. 5 as the measured signal SF which has the magnitude 
of zero centered about v 0 . This data provides information 
about pure Rayleigh scattering that can be used to calculate 
the ratio of Mie scattering to Rayleigh scattering, as shown 
,50 in Eq. 17: 

•S's(i'o) _ r E S D S ,rLFR(v 0 ) + mLF(vo) (Eq. 17) 

Sn(v 0 ) l E n D n \ r + m 


Ss ^ P l T l T r E s D s J dvidv r (L(vi)F(Vi — v))(rR(v — (v r — Av D )) + mM{v — (v r — Aiq>))) 

•Si(v) PiTfT R E\: jd Vj d v r L( Vj )( rR( v - ( v r - Avd)) + mM(v - (v r - Av D ))) 


(Eq. 13) 


Since M is a delta function. Equation 13 simplified to: 

S 5 (v) r EsDs irLFR(v- Av^) + mLF(v - Av^) (Eq. 14) 

Sn(v) r + m 

In one embodiment, tuning laser 141 to a frequency (v A ) 
far enough removed from the effects of the vapor filter 152 
enables the measurement of the ratio of the optical and 
detector efficiencies of the signal channels 158 (Ss as 
represented by Eq. 2 above) and 155 (S N as represented by 
Eq. 5 above). This enables a check for abnormalities in the 
filter on a scan-by-scan basis: 

Ss(v„) _ E s d s (Eq. 15) 

• 5 n ( v „) E n D n 


Since the vapor filter fully attenuates the Mie scattering in 
this region: 

Ss(vq) _ I E s D 5 yLFR(vo) (Eq. 18) 

S N (v 0 ) r + m 

where LFR(v 0 ) is the value of the theoretical return signal at 
particular atmospheric temperature and pressure. 
Accordingly, computer 156 calculates the ratio of Mie 
scattering by first defining a variable Kq as follows: 

„ _ •Ss(vo) (Eq. 19) 

° S«(v 0 ) 

and then solving for the ratio 


45 


50 


55 


In one embodiment, a variable K a can be defined as: 


m 

r 


— LFR{V(y) — 1 


(Eq. 20) 


= Ss(v„) (Eq. 16) 

° S N (v„) 

Once both data sets (i.e., S s and Sv) are symmetric about 
the same data point, computer 156 calculates temperature 


Using the normalized signal return in the region of interest 
65 (i.e., the sloped region between the minimum and maximum 
of the signal return) and writing the result in terms of the 
ratio of m over r, yields the following: 
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S s (v) LFR(v) + ™LF(v) 

S^v)~ Ka 77^_ 


(Eq. 21) 


5 


Substituting the ratio of m and r of Eq. 20 into Eq. 21 yields: 


Sj(v) LFm r Ko i (Eq. 22) 

Sjv(v) Ka LFR(v 0 ) + K a LRF(v 0 )\ 


Solving for LRF(v) yields: 

1 Z.ffF(v 0 n (Eq. 23) 

K„ K 0 J’ 


LFR(v) 


S s (v) LFR(v 0 ) 
Sjv(v) a;. 


EE-(v) 


where the measured signal return LRF(v) is written in terms 
of measured quantities and the theoretical values of LRF 
(v 0 ). Computer 156 then calculates LRF(v) and compares it 20 
to the lookup table 170 to determine atmospheric tempera- 
ture and pressure, described in greater detail in FIG. 5. 

Accounting for power fluctuations, optical efficiencies 
and detector efficiencies, as described herein, allows for an 
independent check on vapor filter 152 during the operation 25 
OADS 140. With variable characteristics of detector chan- 
nels and power fluctuations accounted for, computer 156 
determines, for example, the substantially invariable char- 
acteristics of vapor filter 152, such that more accurate 
measurements of received backscattered laser radiation 30 
(e.g., laser radiation 148) are obtained. 

In one embodiment, the normalization channel depicted in 
FIG. 4 is used to remove atmospheric power fluctuations of 
laser radiation 148. In doing so, computer 156 measures 
Rayleigh and Mie components of laser radiation 147 in 35 
terms of optical efficiencies and detector efficiencies. Such 
efficiencies are typically measured on a shot-by-shot basis 
during the analysis process. In an exemplary operation, laser 
141 generates and transmits laser energy 142 as a series of 
pulses at a particular pulse repetition frequency (“PRF”). 40 
Computer 156 then measures the Rayleigh and Mie com- 
ponents in terms of optical efficiencies and detector effi- 
ciencies on a pulse-by-pulse basis. 

To measure Rayleigh components and Mie components, 
in one embodiment, OADS 140 tunes the frequency of the 45 
laser energy 142 transmitted by laser 141. For example, laser 

141 transmits the laser energy 142 at distal frequencies from 
the peak absorption frequency of filter 152 (illustrated by Vy 
in FIG. 4) to provide a frequency-independent measurement. 
Computer 156 then determines the line shape of laser energy 50 

142 through filter 152. 

In one embodiment, measured intensity of the detected 
backscattered laser radiation (e.g., as determined by elec- 
tronic signal 158) is functionally compared to normalized 
atmospheric factors. Since the measured intensity often 55 
depends upon Mie scatterers (e.g., aerosols) and air density 
changes due to altitude changes and temperature changes. 
The air density changes and the temperature changes are not, 
however, removed through the normalization processes 
described herein; for computer 156 to accurately determine 60 
air parameters such as temperature and pressure of air 144, 
air density changes are removed from the detected back- 
scattered laser radiation so that computer 156 can accurately 
determine the air parameters. 

FIG. 5 shows graph 240 with curves 241 (detected back- 65 
scattered laser radiation at an air density causing Rayleigh 
scattering), 242 (detected backscattered laser radiation at a 
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higher air density causing both Rayleigh and Mie scattering) 
and 243 (theoretically pure Rayleigh scattering). Curves 
241, 242 and 243 illustrate laser energy magnitudes (plotted 
with respect to axis 250) as a function of frequency (plotted 
with respect to axis 251). In one embodiment, computer 156 
processes data from curves 241, 242 and 243 to determine 
other air parameters. For example, Mie scattering effects are 
substantially isolated and removed from calculations to 
determine air temperature and air pressure, since these Mie 
scattering effects produce inaccurate measurements due to 
inconsistent aerosol concentrations. 

To determine the air temperature and air pressure, com- 
puter 156 processes the data from curves 241, 242 and 243 
to substantially isolate and remove the Mie scattering 
effects, such as those found in curve 242. In processing the 
data from curves 241, 242 and 243, computer 156 calculates 
lookup table 170 in substantially real time using a measured 
laser/filter profile (i.e., as measured at detector 165 of FIG. 
2) convolved with theoretical Rayleigh functions for a 
particular temperature and pressure (e.g., illustrated by 
curves 242 and 243). Computer 156 then scales the mea- 
sured return signal LFR(v) (i.e., illustrated by curve 241 in 
this example) with the ratio of m to r determined by Eq. 20. 
Computer 156 then analyzes data near the deepest portion of 
the filter attenuation (i.e., approximately +/-0.5 GHz from 
v f ) to estimate pressure and/or temperature. This portion 
corresponds to a 60 dB region of absorption not contami- 
nated by Mie scattering, a preferred aspect of the invention 
that provides temperature and pressure accuracy by provid- 
ing a reliable temperature base from which to increment 
temperature and/or pressure estimates. 

Computer 156 calculates theoretical Rayleigh return 
assuming an initial temperature estimate and performs a 
Least Square Error (LSE) calculation to determine the 
accuracy of the temperature with respect to the theoretical 
Rayleigh function. Computer 156 repeats the process with 
incremental changes to temperature and/or pressure until an 
optimal fit (i.e., a LSE calculation that corresponds to design 
specifications) is achieved. Although discussed in detail with 
respect to LSE, other approximation methods, such as 
Newton-Rhapson and Monte Carlo, may be used in alter- 
native embodiments. Accordingly the invention is not 
intended to be limited to the calculation methods of the 
embodiments described herein. 

Temperature affects air density in a manner that is recip- 
rocal to pressure; increasing pressure increases density, 
while increasing temperature decreases density. 
Additionally, increasing temperature increases the Rayleigh 
lineshape width while increasing pressure increases the 
Rayleigh lineshape height. Accordingly, for each incremen- 
tal value of temperature and/or pressure, the Rayleigh line- 
shape is unique. Those skilled in the art should readily 
understand that such scattering theory is generally defined in 
“On The Kinetic Model Description Of Rayleigh-Brillouin 
Scattering From Molecular Gases”, G. C. Tenti, D. Boley 
and R. C. Desai, Canadian Journal of Physics, vol. 52, pg. 
285-290 (1974). 

In one example, computer 156 determines air density 
changes by aligning peak absorption frequencies of curves 
241, 242 and 243, illustrated at frequency v f . Since curve 
241 represents detected backscattered laser radiation con- 
taining substantially no Mie scattering, curve 241 can be 
used as a reference where Mie scatterering has been elimi- 
nated. In one example, computer 156, therefore, uses curve 
243 to remove the effects of Mie scattering, aligning curves 
241, 242 and 243 and by calculating a ratio of the detected 
backscattered laser radiation to a theoretically pure Rayleigh 
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scattering which then determines air density (the ratio of 
curves 241 and 242). Mie scattering effects are then removed 
by subtracting curve 243 from the calculated ratio of curves 
241 and 242. With the Mie scattering essentially removed 
from the measurement, computer 156 more accurately deter- 
mines air temperatures and air pressures. 

FIGS. 6 and 7 show other exemplary graphs that may be 
used in determining air pressure and air temperature. FIG. 6 
illustrates a graph 260 of electronic signals 163 and 167 of 
FIG. 2 respectively generated by detectors 162 and 165 of 
FIG. 2. Graph 260 shows the electronic signals 163/167 
representing light intensity as a function of normalized 
signal strength (axis 261) versus frequency (axis 262). FIG. 
7 illustrates a graph 280 of electronic signals 158 and 155 of 
FIG. 2 respectively generated by detectors 153 and 154 and 
representing light intensity as a function of normalized 
signal strength (axis 281), versus frequency (axis 282). 
These four light intensities (represented by electronic signals 
163, 167, 158 and 155) are measured, over time, through a 
collection of laser pulses. For example, the laser energy 142 
of FIG. 2 generated by laser 141 at a certain PRF is swept 
in transmitted frequency such that each laser pulse is emitted 
at a different frequency. Electronic signals 163 and 167 
therefore illustrate how laser 141 may sweep laser energy 
142 in frequency across an absorption band 263 of the vapor 
filter 152. Illustratively, FIG. 6 shows one complete fre- 
quency sweep of laser energy 142 generated by laser 141 
and detected by detectors 162 and 165. Similarly, electronic 
signals 155 and 158 of FIG. 7 show detected signals of 
detectors 162 and 165 (backscattered) as laser 141 performs 
a complete sweep in frequency of laser energy 142 across 
the absorption band 283 of the vapor filter 152. 

From signals 163 and 167, computer 156 may for example 
determine a normalized filter transmission, by dividing 
discrete points of electronic signal 167 by corresponding 
discrete points of signal 163. Similarly, computer 156 may 
determine a normalized atmospheric return though vapor 
filter 152 by dividing discrete points of signal 158 by 
corresponding discrete points of signal 155. These discrete 
points, described herein, correspond to individual pulses of 
laser energy 142. 

Using normalized calculations of filter transmission (e.g., 
from graph 260) and the normalized calculations of atmo- 
spheric return (e.g., from graph 280), computer 156 deter- 
mines relative optical efficiencies in the vapor filter 152. 

In one embodiment, computer 156 determines optical 
transmission for vapor filter 152 using the frequency inde- 
pendent components of data from graph 260; FIG. 6 (there 
is substantially no change in amplitude for signals 163 and 
167 at frequencies greater in magnitude than ±18 GHz from 
0 GHz illustrated at points 264, 265, 266 and 267). Com- 
puter 156 therefore determines the ratio of optical transmis- 
sion for vapor filter 152 by calculating a ratio of signal 167 
to signal 163, via frequency corresponding points of the 
signals, for points representing frequencies greater in mag- 
nitude than ±18 GHz from 0 GHz. 

Similarly, computer 156 determines a magnitude of inten- 
sity of atmospheric-returned laser energy received through 
vapor filter 152 using the frequency independent parts of the 
data from graph 280 of FIG. 7 (there is substantially no 
change in amplitude for signals 155 and 158 at frequencies 
greater in magnitude than ±18 GHz from 0 GHz illustrated 
at points 284, 285, 286 and 287). Computer 156 thereby 
determines a ratio of atmospheric return with the laser power 
measurement by calculating a ratio of signal 158 to signal 
155 via frequency corresponding points of the signals for 
points representing the frequencies greater in magnitude 
than ±18 GHz from 0 GHz. 
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In one embodiment, computer 156 calculates a ratio of 
signal 158 to signal 155 for frequencies between ±0.5 GHz 
(illustrated at points 288 and 289). Such a frequency range 
substantially no Mie scattering of laser energy 142 for air 
s 144; it thus indicates substantially pure Rayleigh scattering. 
Computer 156 thus compares a Rayleigh to Mie scattering 
strength based upon the ratio of signal 158 to signal 155. 
Computer 156 determines Rayleigh to Mie scattering 
strength by comparing the ratio of signal 158 to signal 155 
at frequencies between ±0.5 GHz to the ratio of signal 158 
to signal 155 at frequencies greater than ±18 GHz from 0 
GHz. In one embodiment, computer 156 performs similar 
calculations for "non-scattered” laser energy 142 (e.g., com- 
ponent 143B of FIG. 2) based on data illustrated in FIG. 6 
using points 268 and 269. Such a process is further described 
15 in FIG. 11. 

The ratios determined for the non-scattered laser energy 
142 and for the scattered laser energy 142 may be used in 
tandem to numerically calculate Rayleigh-Laser-Filter con- 
volution from data. The Rayleigh-Laser-Filter convolution 
20 is in turn compared to a look up table of theoretical 
Rayleigh-Laser-Filter convolution values to determine tem- 
perature and pressure. 

FIGS. 8 and 9 represent errors in pressure and temperature 
measurements obtained during experimentation. FIG. 6 
25 shows the percentage error in pressure measurements as a 
function of pressure and temperature, obtained using syn- 
thetic data with a 0.5% noise component. FIG. 9 shows the 
percentage error in temperature measurements as a function 
of pressure and temperature obtained using synthetic data 
.to with a 0.5% noise component. 

FIG. 10 shows one OADS transceiver 300. OADS trans- 
ceiver 300 is for example suitable for use as transceiver 110 
of FIG. 1. OADS transceiver 300 is illustrated in FIG. 10 as 
a midsectional view to reveal interior components of the 
35 transceiver. In one embodiment, OADS transceiver 300 is 
hermetically sealed and configured for projecting laser 
radiation 301 to air 302 and for detecting backscattered laser 
radiation 303. Optical coupling 305 connects to optical fiber 
306. Optical fiber 306 communicates laser energy radiated 
40 from a laser to project laser radiation 301 to air 302. OADS 
transceiver 300 may include beam steering optics 307 and 
lens 308 to assist in the projection of radiation 301. 

OADS transceiver 300 receives backscattered laser radia- 
tion 303 when laser radiation 301 backscatters from air 302. 
45 Transceiver 300 receives backscattered laser radiation 303 
through a receiving aperture 309 embedded with barrier 310. 
Aperture 309 captures and directs backscattered laser radia- 
tion 303 to rigidly mounted mirrors 311 (mounted within 
transceiver 300 by rigid mounting brackets 320) for 
50 improved detection. Mirrors 311 may be annular focusing 
mirrors configured for demagnifying the backscattered laser 
radiation 303 received through aperture 309. The rigid 
mounting brackets substantially prevent jitter of backscat- 
tered laser radiation 303. 

55 Demagnified laser radiation 312 is reflected towards fold- 
ing mirrors 313, which, in turn, direct laser radiation 312 to 
partial reflector 314. Partial reflector 314 reflects a portion of 
laser radiation 312 in the form of laser radiation 315 towards 
detector 316. Partial reflector 314 is a beam splitter that 
60 reflects a portion (e.g., 50%) of laser radiation 312 and 
transmits a portion (e.g., 50%) of laser radiation 312. 

Detector 316 is similar to detector 154 of FIG. 2. As such, 
detector 316 converts received laser radiation 315 into an 
electronic signal, such as electronic signal 155 of FIG. 2, for 
65 computer processing. A computer, such as computer 156 of 
FIG. 2 processes the electronic signal to remove laser power 
fluctuations caused by atmospheric changes in air 302. 
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Laser radiation 312 that is directed to vapor filter 317 is 
filtered by vapor filter 317. In one embodiment, vapor filter 
317 is an atomic vapor filter or a mercury vapor filter and 
electronically controlled by electronics 318. Electronics 318 
may control temperature and pressure of the vapor within 
filter 317 to stabilize or electronically “lock” frequency 
characteristics of filter 317. Such stabilization makes OADS 
transceiver 300 less susceptible to vibrational misalignment. 

Laser radiation 312 is filtered by filter 317 and passes to 
detector 319 in the form of filtered laser radiation 321. 
Detector 319 receives and converts the filtered laser radia- 
tion 321 into an electronic signal, which is also processed by 
the computer to determine the desired air parameters. 

FIG. 11 shows a flowchart of one exemplary methodical 
operation 400 of an OADS. A laser (e.g., laser 141 of FIG. 
2) sweeps laser energy across a predetermined frequency 
spectrum, in element 401. The laser may sweep the laser 
energy across a frequency range of about +/-20 GHz by 
firing the laser energy at a certain PRF. In one embodiment, 
the PRF is about 1 kHz, with a pulse width between about 
50 ns and 100 ns, and the swept frequency range is centered 
about a frequency corresponding to a peak absorption fre- 
quency (e.g., 260 nm) of a vapor filter (e.g., vapor filter 152, 
FIG. 2). 

Laser energy is typically split into a four distinct paths 
such that the laser energy can be detected as four different 
inputs, in element 402. These four paths of laser energy 
correspond to: 1) laser energy transmitted by the laser (e.g., 
component 159 of FIG. 2); 2) laser energy transmitted by the 
laser through the vapor filter (e.g., component 164 of FIG. 
2); 3) laser energy transmitted by the laser into the air and 
backscattered (e.g., component 150Aof FIG. 2); and 4) laser 
energy transmitted by the laser into the air and backscattered 
through the vapor filter (e.g., component 157 of FIG. 2). For 
simplicity, these components are hereinafter referred to as: 
1) unfiltered laser energy; 2) filtered laser energy; 3) unfil- 
tered backscattered laser energy; and 4) filtered backscat- 
tered laser energy. 

After detecting the four forms of laser energy, a computer 
(e.g., computer 156, FIG. 2), determines normalized filter 
transmission of the vapor filter, in element 403. For example, 
the computer in one embodiment processes the unfiltered 
laser energy and filtered laser energy by dividing the mag- 
nitude of the filtered laser energy by the magnitude of the 
unfiltered laser energy. The division is performed on a pulse 
by pulse basis, where divided magnitudes of the pulses have 
corresponding frequencies. 

The computer also determines, in one embodiment, nor- 
malized atmospheric return of the laser energy, in element 
404. For example, the computer processes the unfiltered 
backscattered laser energy and unfiltered backscattered laser 
energy by dividing the magnitude of the filtered backscat- 
tered laser energy by the magnitude of the unfiltered back- 
scattered laser energy. Again, division is performed on a 
pulse by pulse basis, where divided magnitudes of the pulses 
have corresponding frequencies. 

Once normalized filter transmission and normalized 
atmospheric return of the laser energy are determined, the 
computer determines signal strengths for each of the filter 
transmission and the atmospheric return. For example, the 
computer determines the optical transmission through the 
filter by calculating a ratio of the filtered laser energy to the 
unfiltered laser energy at particular frequency ranges, in 
elements 405 and 407. The computer similarly determines 
the atmospheric return (scattering) signal strength through 
the filter by calculating a ratio of the filtered backscattered 
laser energy to the unfiltered laser energy at particular 
frequency ranges, in elements 406 and 408. 
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The computer also determines a signal strength ratio for 
the normalized filter transmission by dividing filtered laser 
energy by unfiltered laser energy, again on a pulse by pulse 
basis, at frequencies greater in magnitude than about +/-18 
5 GHz about the peak absorption frequency, in element 407. 
The computer further determines a signal strength ratio for 
the normalized filter transmission by dividing filtered laser 
energy by unfiltered laser energy on a pulse by pulse basis 
at frequencies between about +/-0.5 GHz, in element 405. 
10 These signal strength determinations correspond to fre- 
quency ranges where Mie scattering (e.g., +/-18 GHz) and 
Rayleigh scattering (e.g., +/-0.5 GHz) are most prevalent — 
useful when combined with similar signal strength determi- 
|S nations for the normalized atmospheric return. The com- 
puter determines a Mie scattering signal strength ratio for 
the normalized atmospheric return of the laser energy by 
dividing filtered backscattered laser energy by unfiltered 
backscattered laser energy, again on a pulse by pulse basis, 
20 at frequencies greater in magnitude than about +/-18 GHz 
about the peak absorption frequency, in element 408. The 
computer also determines a Rayleigh scattering signal 
strength ratio for the normalized atmospheric return of the 
laser energy by dividing filtered scattered laser energy by 
25 unfiltered backscattered laser energy on a pulse by pulse 
basis at frequencies between about +/-0.5 GHz in element 
406. 

With signal optical transmission for the filter and signal 
strengths for both Rayleigh scattering and Mie scattering 
determined, the computer determines a Rayleigh laser filter 
convolution, in element 409. For example, the computer in 
one embodiment performs a convolution of the optical 
transmission with the Rayleigh and Mie scattering signal 
, 5 strengths corresponding to the frequency ranges for Ray- 
leigh and Mie scattering of +/-0.5 GHz and +/-18 GHz, 
respectively. The computer then accesses a lookup table, 
such as lookup table 170 of FIG. 2, having theoretical 
Rayleigh laser filter convolution values to determine tem- 
4(] perature and pressure of the air, in element 410. 

Advantages to the embodiments described above include 
improved computations of various air parameters, such as 
air speed, air temperature and air pressure, substantially 
regardless of altitude and/or Mie scattering. Other advan- 
45 tages include a system that accurately performs in a variety 
of vibrational environments. Still, other advantages include 
an ability to determine temperature and pressure within a 
particular region of atmosphere without a prior knowledge 
of the atmosphere. Additionally, the instrumentalities as 
50 shown and described above advantageously require reduced 
on-aircraft system calibrations, a more robust system in 
vibrational environments, and system health checks. 

Since certain changes may be made in the above methods 
and systems without departing from the scope of the 
55 invention, one intention is that all matter contained in the 
above description or shown in the accompanying drawings 
be interpreted as illustrative and not in a limiting sense. By 
way of example, those skilled in the art should appreciate 
that the OADS and the OADS transceivers, as described 
60 herein, may be constructed, connected, arranged, and/or 
combined in other various manners without departing from 
the scope of the invention. Another intention includes an 
understanding that the following claims are to cover all 
generic and specific features of the invention described 
65 herein, and all statements of the scope of the invention 
which, as a matter of language, might be said to fall there 
between. 
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What is claimed is: 

1. A method for remotely sensing air outside a moving 
aircraft, comprising: 

projecting laser radiation from the aircraft into the air to 
induce scattered radiation that has a molecular scattered 5 
radiation component and an aerosol scattered radiation 
component; 

detecting scattered laser radiation; 

distinguishing the molecular scattered laser radiation 10 
component from the aerosol scattered radiation com- 
ponent; and 

determining one or more air parameters based on the 
scattered radiation, the air parameters selected from the 
group of air speed of the aircraft, pressure outside the is 
aircraft, temperature outside the aircraft and aircraft 
orientation angle. 

2. A method of claim 1, wherein distinguishing comprises 

automatically analyzing a spectrum of the molecular and 
aerosol scattered radiation components. 20 

3. A method of claim 2, wherein analyzing comprises 
measuring a spectral lineshape of the spectrum. 

4. A method of claim 1, wherein distinguishing comprises 
distinguishing Rayleigh characteristics of the scattered 
radiation from Mie characteristics of the scattered laser 25 
radiation. 

5. A method of claim 1, the step of determining one or 
more air parameters comprising determining air speed, 
wherein determining air speed comprises determining a 
Doppler line shift of a spectral line shape of at least one of .50 
the molecular scattered radiation and the aerosol scattered 
radiation. 

6. A method of claim 1, the step of determining one or 
more air parameters comprising determining air speed and 
aircraft orientation angle, and further comprising determin- 35 
ing the air speed and the aircraft orientation angle using a 
molecular scattered Doppler line shift and without using 
aerosol distribution. 

7. A method of claim 5, wherein an aerosol scattered 
Doppler line shift is used in determining the air speed and 40 
the aircraft orientation angle. 

8 . A method of claim 1, the step of determining one or 
more air parameters comprising determining temperature, 
wherein determining temperature comprises determining a 
spectral line shape of the molecular scattered laser radiation 45 
and comparing the line shape to one or more theoretical 
Rayleigh line shapes. 

9. A method of claim 1, the step of determining one or 
more air parameters comprising determining pressure, 
wherein determining pressure comprises determining ampli- 50 
tude of a line shape from the molecular scattered radiation. 

10. A method of claim 1, the step of determining one or 
more air parameters comprising determining pressure, 
wherein determining pressure comprises determining a line 
shape from the molecular scattered radiation and comparing 55 
the line shape to one or more theoretical Rayleigh line 
shapes. 

11. A method of claim 1, wherein projecting laser radia- 

tion comprises utilizing one or more tunable narrow lin- 
ewidth lasers. 60 

12. A method of claim 11, the lasers having a center 
frequency matched to one of a mercury vapor filter peak 
absorption frequency and an atomic vapor filter peak absorp- 
tion frequency. 

13. A method of claim 1, wherein projecting comprises 65 
dividing the laser radiation through a plurality of transceiv- 
ers and optical fibers. 


18 

14. A method of claim 13, further comprising utilizing 
two or more transceivers mounted with the aircraft, wherein 
detecting comprises detecting the scattered radiation with 
the two or more transceivers. 

15. A method of claim 1, wherein detecting scattered laser 
radiation comprises utilizing at least one of a fixed fre- 
quency atomic vapor filter and a fixed frequency mercury 
vapor filter. 

16. A method of claim 1, wherein distinguishing com- 
prises deconvolving Rayleigh line shapes and Mie line 
shapes via digital signal processing. 

17. A method of claim 1, wherein a spectral lineshape of 
the molecular scattered radiation component is used to 
determine the temperature and the pressure. 

18. A system for sensing of air outside a moving aircraft, 
comprising: 

at least one laser for generating laser energy; 

at least one transceiver for projecting the laser energy 
from the aircraft to the air and for receiving scattered 
laser energy from the air; and 

a computer for processing signals from the transceiver to 
distinguish molecular scattered radiation from aerosol 
scattered radiation and for determining one or more air 
parameters based on the scattered laser radiation. 

19. A system of claim 18, the parameters selected from a 
group consisting of air speed of the aircraft, pressure outside 
the aircraft, temperature outside the aircraft and aircraft 
orientation angles. 

20. A system of claim 18, further comprising at least one 
fixed frequency atomic vapor filter. 

21. A system of claim 20, the laser comprising a tunable 
frequency centered at an absorption wavelength of said filter. 

22. A system of claim 18, the at least one atomic vapor 
filter comprising one of mercury or cesium. 

23. A system of claim 18, the laser energy comprising a 
wavelength in a range of about 254 nm to 355 nm. 

24. A system of claim 18, the computer comprising means 
for deconvolving Rayleigh line shapes and Mie line shapes 
of the molecular and the aerosol scattered radiation. 

25. A system of claim 18, further comprising optical fiber 
for coupling the laser to the transceiver. 

26. A system of claim 18, further comprising a transducer 
for converting scattered laser radiation to electronic signals 
for the computer. 

27. A system of claim 18, the system further comprising 
means for determining a temperature and a pressure within 
a region of the air without a prior knowledge of the air. 

28. A transceiver, comprising: 

beam steering optics for projecting laser energy to air; and 

a vapor filter configured for filtering backscattered laser 
energy received from the air, wherein filtered backscat- 
tered laser energy is processable to determine two or 
more air parameters selected from a group consisting of 
air speed, pressure and temperature. 

29. A transceiver of claim 28, further comprising at least 
one detector configured for receiving the filtered backscat- 
tered laser energy and for converting the filtered backscat- 
tered laser energy to electronic signals for use in determining 
the one or more air parameters. 

30. A transceiver of claim 28, comprising a plurality of 
mirrors configured for directing the laser energy through the 
vapor filter. 

31. A transceiver of claim 28, the beam steering optics 
comprising an optical connector configured for coupling to 
a laser generating the laser energy. 
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32. A transceiver of claim 28, comprising one or more 
mounts configured for mounting the transceiver within an 
aircraft. 

33. A transceiver of claim 28, the beam steering optics 
comprising a lens. 

34. A transceiver of claim 28, wherein filtered backscat- 
tered laser energy is processable to determine aircraft ori- 
entation angles. 
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35. A transceiver of claim 28, the vapor filter comprising 
a fixed frequency atomic vapor filter or a fixed frequency 
mercury vapor filter. 

36. A transceiver of claim 28, the laser energy comprising 
5 a center frequency corresponding to a peak absorption 

frequency of the vapor filter. 
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